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Abstract 

A  kinetic  model  of  a  laser-induced  breakdown  spectroscopy  (LIBS)  plume  of  cyclotrimethylenetrinitramine  (RDX)  was  developed  for  the 
analysis  of  processes  responsible  for  the  LIBS  signature  of  explosives.  Air  and  argon  were  considered  as  buffer  gases.  The  model  includes  a  set  of 
processes  involving  ion  chemistry,  as  well  as  excitation,  ionization,  and  other  processes  affecting  neutral  and  ion  species.  Modeling  results  show 
that  the  overall  reaction  process  can  be  considered  a  two-stage  process.  The  first  stage  corresponds  to  a  fast  approach  to  a  quasi-stationary  state, 
while  the  second  stage  corresponds  to  the  change  of  quasi-stationary  species  concentrations  due  to  the  change  in  temperature.  As  a  result  of  the 
two-stage  process,  the  initial  mechanism  of  explosive  decomposition  is  not  important  in  determining  its  signature  in  the  LIBS  measurement  time 
window  (1-30  ps).  The  main  processes  responsible  for  generation  of  excited  states  for  the  LIBS  emission  are  electron-excitation  impact  processes. 
A  mechanism  for  the  appearance  of  a  double  peak  of  the  C2  species  concentration  in  the  RDX  plasma  plume  was  suggested.  Double-peak  behavior 
of  the  C2  species  was  previously  experimentally  observed  during  laser  ablation  of  graphite. 
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1.  Introduction 

Laser-induced  breakdown  spectroscopy  (LIBS)  is  an  emer¬ 
ging  diagnostic  method  for  many  applications.  A  LIBS  operates 
with  a  high  power  laser  pulse  focused  on  small  spot  of  a  sample 
material.  The  interaction  of  the  pulsed  laser  beam  with  a  sample 
material  produces  high-temperature  ionized  plasma  containing 
electronically  excited  elements  which  radiate  the  characteristic 
emission  lines  of  corresponding  elements.  It  is  important  to  note 
that  the  laser  energy  release  time  occurs  in  the  femtosecond- 
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nanosecond  temporal  range.  However  a  plasma  kernel  with  high 
temperature  exists  for  10-100  ps,  as  a  result  of  chemical 
reactions  and  recombination  processes,  releasing  energy  on  the 
microsecond  time  scale  and  sustaining  the  plasma. 

The  direct  detection  of  energetic  materials  and  explosives  in 
real  time  is  an  important  practical  problem  [1-3].  The  LIBS 
technique  is  mainly  considered  to  be  a  method  of  analytical 
element  analysis.  However,  the  LIBS  technique  allows  one  to 
identify  certain  chemical  compounds  with  chemometric  analysis. 
Results  of  recent  LIBS  studies  of  energetic  materials  demonstrate 
that  the  simultaneous  registration  O,  N,  H  and  C2  signals  [1-3] 
and  analysis  of  the  ratios  of  peak  intensities  allows  one  to  make 
definite  conclusions  on  the  presence  of  energetic  material. 

The  aim  of  the  present  work  was  to  develop  a  kinetic  model  of 
the  LIBS  plume  of  RDX  for  the  analysis  of  processes  responsible 
for  the  RDX  LIBS  signature.  RDX  (1,3,5-trinitrohexahydro-s- 
triazine)  is  an  explosive  nitroamine  widely  used  in  military  and 
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industrial  applications.  It  is  considered  one  of  the  most  powerful 
military  explosives.  The  developed  model  includes  processes  of 
fonnation  and  quenching  of  excited  H,  N  and  O  atoms,  and  a 
submodel  for  formation  of  the  C2  molecule.  The  observation  of 
LIBS  C2  signals  was  used  to  discriminate  organic  samples  from 
inorganic  targets  during  detection  of  energetic  materials  [2], 
Together  with  peak  ratios,  the  data  on  C2  emission  was  used  for 
the  detection  of  a  number  of  explosive  materials.  In  future  studies, 
we  plan  to  reduce  suitably  the  developed  model  for  detailed  CFD 
modeling  studies  of  RDX  LIBS  plume.  The  results  of  the  present 
kinetic  modeling  suggest  that  the  assembled  model  can  be  used  for 
understanding  the  reaction  behavior  in  the  LIBS  plume  of  other  C / 
H/N/O-containing  explosive  materials. 

2.  Kinetic  model 

As  the  first  step  in  the  analysis  of  processes  responsible  for 
the  emission  lines  in  the  plasma  plume  of  explosives,  a  kinetic 
model  of  RDX  reaction  behavior  in  the  LIBS  plume  was 
developed.  Air  and  argon  were  considered  as  buffer  gases. 
Earlier  we  developed  a  kinetic  model  of  lead  reactions  in  air  and 
argon  [4,5]  to  describe  the  chemical  reactions  and  collisional 
processes  that  occur  within  the  high-temperature  mixture  in  the 
LIBS  plume.  These  processes  were  included  as  a  part  of  the 
current  model.  The  model  includes  a  set  of  processes  involving 
ion  chemistry,  as  well  as  excitation,  ionization,  and  other 
processes  affecting  neutral  and  ion  species. 

The  kinetic  model  of  RDX  reactions  in  the  plume  includes 
the  following  main  reaction  subsets:  RDX  decomposition 
reactions,  high-temperature  air  and  argon  plasma  reactions  of 
neutral  species,  reactions  of  ionization  and  reactions  of  charged 


Table  1 


States  of  atomic  nitrogen,  hydrogen  and  oxygen  included  in  the  kinetic  model 


Atom 

Level 

Energy,  cm *  1 

Emission  line,  mn 

N 

2p3  4S° 

0 

(L)520, 347, 120, 113.4 

N*2 

2p3  2D° 

19,226 

(U)520,(L)149 

N*3 

2p3  2P° 

28,839 

(U)347,(L)175 

N*4 

2p23s 4P 

83,337 

(U)  1 20,(L)(870),(82 1 ), 747,742 

N*5 

2p23s 2P 

86,192 

(U)149,175,(L)135,947,(940),742.861 

N*6 

2s22p4  4P 

88,134 

(U)  113 

N*7 

2p23p  2S° 

93,582 

(L)  947,939,(857-866) 

N*8 

2p23p  4D° 

94,839 

(U)(868-874) 

N*9 

2p23p  4P° 

95,511 

(U)821 

N*10 

2p23p  4S° 

96,752 

(U)745,947 

N*ll 

2p23p  2D° 

96,834 

(U)740.9, 939.5 

N*12 

2p23p  2P° 

97,793 

(U)861 

0 

2p4  3P 

0 

0*2 

2p4  ‘D 

15,865 

0*3 

2p4  ‘S 

33,795 

0*4 

3s  5S° 

73,769 

(L)777 

0*5 

3s  3S° 

76,794 

(L)844 

0*6 

3p  5P 

86,626 

(U)777 

0*7 

3p  3P 

88,642 

(U)844 

0*8 

4s  5S° 

95,498 

H 

n=l 

0 

(L)102.5,  121.6 

H*2 

n=2 

82,259 

(L)656,(U)102.5 

H*3 

n=3 

97,492 

(U)121.6,  656 

Table  2 

Kinetic  models  of  RDX  decomposition  used  in  this  work 
Kinetic  model  Model  details 


Gas-phase  decomposition  model 

Initial  gas  composition  on  the 
surface  from  RDX  decomposition 
and  detailed  gas-phase  model  of 
this  work 

Two  overall  subsurface  reactions  and 
gas-phase  mechanism  of  this  work 

Overall  reaction  of  decomposition  and 
detailed  mechanism  of  this  work 


Model  of  initial  RDX  decomposition 
of  Chakraborthy  et  al.  [8] 
Decomposition  products  at  the  surface, 
experimental  data  of  Fetherolf  and 
Litzinger  [11] 

Subsurface  decomposition  model  from 
[12] 

Overall  decomposition  reaction  from 
[10] 


species  (e-,  H+,  0+,  Oj,  N+,  NO+,  Nj,  Hj,  OH+,  H20+,  O2 , 
O',  C+,  CO+,  CO?,  CH3,  CH?,  CH+,  CH?,  C2H?,  C2H2,  N?, 
N20?,  N30+,  N20+,  N?,  NO?,  NO?,  NO',  N20',  N02 ,  NO?), 
reactions  generating  the  electronically  excited  atoms  N,  O,  and 
H  and  their  reactions  in  the  plasma  plume.  Table  1  contains  the 
electronic  states  of  atomic  N,  O,  and  H  included  in  the  model. 
For  example,  the  following  excited  states  of  the  oxygen  atom 
were  included  in  the  kinetic  model:  ground  state  0(3P),  two 
metastable  states  0(*D)  and  0( 1 S)  and  excited  states  3s  5S,  3s 
3S,  3p  5P,  and  3p  3P,  which  participate  in  the  777  nm  and 
844  nm  O  atom  transitions. 

Thus  the  model  includes  the  following  reaction  subsets: 

1.  Initial  decomposition  reactions  of  RDX  (or  explosive 
products  at  the  surface  of  target  material) 

2.  High-temperature  reactions  of  products  of  explosive  decom¬ 
position  in  the  buffer  gas  under  the  LIBS  conditions  (buffer 
gas:  Ar  or  air). 

3.  Reactions  of  excited  states  species  responsible  for  the  LIBS 
signature,  including  the  generation  and  quenching  of  excited 
states  (O,  H,  N)  and  the  subset  of  reactions  of  excited  species 
with  the  decomposition  products  of  RDX  and  with  buffer  gas 
species. 

The  kinetic  model  was  assembled  from  several  sources.  Main 
sources  of  data  are  identified  below.  When  absent,  kinetic  data 
were  estimated  based  on  analogy  or  with  the  use  of  empirical 
rules.  The  kinetic  model  contains  overall  137  species,  with  577 
reactions.  To  describe  RDX  decomposition,  several  mechanisms 
were  used  [6- 1 0]  (Table  2).  A  plasma  kinetic  model  for  the  H/N/ 
O/Ar  system  was  taken  from  our  previous  work  [4,5],  Plasma 
chemistry  reactions  related  to  the  C/H/O/N  system  came  from 
several  sources  [13-15].  Reactions  involving  the  C2  species 
were  added  on  the  basis  of  the  results  of  Kruse  and  Ross  [  1 6, 1 7] . 
Kinetic  data  for  electronically  excited  O,  H  and  N  atoms 
(electron-impact  excitation  cross  sections,  radiative  quenching 
rates,  etc)  were  taken  from  [18-21],  Table  3  contains  the  block  of 
reactions  describing  the  kinetics  of  C2  formation  and  consump¬ 
tion,  and  a  block  of  ion  reactions  involving  carbon-containing 
species.  The  kinetic  submodel  describing  the  formation  and 
decay  of  electronically  excited  H,  N  and  O  atoms  will  be 
presented  in  a  separate  publication.  Note  that  the  kinetic  model 
doesn’t  currently  include  excited  C  atom  states.  The  excited 
states  of  carbon  atom  will  be  considered  in  a  future  work. 

Of  course,  at  these  high  temperatures,  a  large  number  of  rate 
constant  values  was  obtained  by  extrapolation  from  rate 
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Table  3 

Reactions  of  C2  formation  and  carbon-containing  ions  [the  rate  constant  k  (cm, 
s,  mol  units)  is  expressed  in  the  Chemkin  format  as  AT"exp(— Is/RT)] 


No  Reaction 

y4(mole,  cm,  s) 

n 

£(cal/mole)  Reference 

Reactions  of  C2  formation 

1.  c2  +  m=c+c  +  m 

1.50E+ 16 

0.0 

142,400.0 

[17] 

2.  C2+H2=C2H+H 

6.60E+ 13 

0.0 

8000.0 

[17] 

3.  C2H  +  M=C2  +  H+M 

1.74E+35 

-5.16 

114,000.0 

[17] 

4.  C2H  +  C2H=C2H2+C2 

1.80E+ 12 

0.0 

0.0 

[22] 

5.  CN+C=C2+N 

3.00E+ 14 

0.0 

35,968.0 

[23] 

6.  C2+N2=CN+CN 

6.40E+ 1 1 

0.0 

24,900.0 

[16] 

7.  C2+N2=C2N+N 

3.30E+ 14 

0.0 

67,660.0 

[16] 

8.  C2+C2=C3+C 

3.20E+ 14 

0.0 

0.0 

[17] 

9.  H+CH=C+H, 

1.10E+14 

0.0 

0.0 

[8] 

10.  C+C2H=C3  +  H 

7.00E+ 14 

0.0 

0.0 

[17] 

11.  C2  +  CH=C3  +  H 

2.00E+ 14 

0.0 

0.0 

[17] 

12.  CH+CH=C2+H  +  H 

1.00E+14 

0.0 

0.0 

[17] 

13.  CH+CH=C2H2 

1.20E+ 14 

0.0 

0.0 

[14] 

14.  C3H+H=C3  +  H2 

1.00E+14 

0.0 

0.0 

[17] 

15.  C3H2+H=C3H+H2 

6.31E+13 

0.0 

0.0 

[17] 

16.  C3+M=C2+C+M 

4.00E  + 16 

0.0 

150,000.0 

[17] 

17.  CH+M=C+H  +  M 

1.90E+ 14 

0.0 

66,968.0 

[17] 

18.  CH2+M=C+H2+M 

1.15E+ 14 

0.0 

55,820.0 

[17] 

19.  CH2+M=CH+H+M 

2.88E+ 14 

0.0 

79,030.0 

[17] 

20.  C3H+M=C3  +  H+M 

1.58E+ 14 

0.0 

39,420.0 

[17] 

21.  C3H2  +  M=C3H  +  H+M 

1.00E+15 

0.0 

60,000.0 

[17] 

22.  C2  +  0=CO+C 

7.60E+ 14 

0.0 

12,380.0 

[17] 

23.  NO  +  C2=C2N  +  0 

4.70E+ 13 

0.0 

8644.0 

[16] 

24.  NO  +  C2=C20+N 

2.80E+ 13 

0.0 

8644.0 

[16] 

25.  C3  +  0=C0+C2 

3.00E+ 14 

0.0 

0.0 

[16] 

26.  C3  +  0=C20  +  C 

1.00E+14 

0.0 

0.0 

[16] 

27.  C+NO=CN+0 

1.90E+ 13 

0.0 

0.0 

[8] 

28.  C+NO=CO+N 

2.90E+ 13 

0.0 

0.0 

[8] 

29.  C+N2=CN+N 

6.30E+ 13 

0.0 

46,020.0 

[8] 

30.  CN+NO=NCO+N 

9.64E+ 13 

0.0 

42,130.0 

[16] 

31.  CN+NO=N2  +  CO 

2.51E+ 10 

0.0 

0.0 

[16] 

32.  CN+M=C+N  +  M 

2.50E+ 14 

0.0 

141,010.0 

[16] 

33.  C2N+N=CN+CN 

6.03E+ 13 

0.0 

0.0 

[16] 

34.  C2N+0=C20+N 

1.00E+14 

0.0 

0.0 

[16] 

35.  C2N+M=CN+C+M 

1.00E+16 

0.0 

130,000.0 

[16] 

36.  C  +  CO+M=C20+M 

2.29E  + 16 

0.0 

0.0 

[16] 

37.  CN+0=CO+N 

7.70E+ 13 

0.0 

0.0 

[8] 

38.  CN+H2=HCN+H 

2.10E+ 13 

0.0 

4710.00 

[8] 

39.  02+C2=C0  +  C0 

6.60E+ 12 

0.0 

757.0 

[24] 

40.  OH  +  C2=CO  +  CH 

5.00E+ 12 

0.0 

0.0 

[25] 

41.  CN+C02=NCO+CO 

4.00E+ 14 

0.0 

38,400.0 

[26] 

42.  CN+CO=NCO  +  C 

1.50E  + 16 

-0.487 

131,600.0 

[26] 

43.  C0  +  0=>C02 

2.50E+06 

0.0 

3200.0 

[15] 

(radiative  recombination) 

44.  C02+C=CO+CO 

6.00E+08 

0.0 

0.0 

[15] 

45.  C  +  0+M=CO  +  M 

7.26E+ 13 

0.0 

0.0 

[27] 

46.  OH  +  C=H  +  CO 

5.00E+ 13 

0.0 

0.0 

[8] 

47.  C  +  02=0+C0 

5.80E+ 13 

0.0 

576.0 

[8] 

48.  C+CH2=H+C2H 

5.00E+ 13 

0.0 

0.0 

[8] 

49.  C+CH3=H+C2H2 

5.00E+ 13 

0.0 

0.0 

[8] 

50.  C+N2=CN+N 

6.30E+ 13 

0.0 

46,020.0 

[8] 

Reactions  of  C-containing  ions 

1.  C  +  0=E  +  CO+  8.80E+08 

1.0 

66,200.0 

[26] 

2.  NO++C=NO+C+ 

1.00E+13 

0.0 

46,400.0 

[26] 

3.  CO+C+=C+CO+ 

1.00E+13 

0.0 

62,800.0 

[26] 

4.  02+C+=C  +  Oj 

1.00E+13 

0.0 

18,800.0 

[26] 

5.  C+E=E+E+C+ 

3.90E+33 

-3.78 

261,400.0 

[26] 

6.  E+O+^O 

1.07E+ 1 1 

-0.52 

0.0 

[26] 

(radiative  recombination) 

7.  E+C+=>C 

2.02E+ 1 1 

-0.46 

0.0 

[26] 

(radiative  recombination) 

Table  3  ( continued ) 


No  Reaction 

y4(mole,  cm,  s) 

n 

Efcal/mole)  Reference 

Reactions  of  C-containing  ions 

8.  E+CHj=CH+H  +  H  4.30E+ 14 

0.5 

0.0 

[14] 

9.  0+NO+=N+Oj 

7.20E+12 

0.29 

97,200.0 

[26] 

10.  CH4+CHj=CH3  +  CHj 

8.20E+13 

0.0 

0.0 

[28] 

11.  HJ  CHi=H+CHi 

9.60E+ 14 

0.0 

0.0 

[28] 

12.  H2+CtT=H+CHi 

7.20E+ 14 

0.0 

0.0 

[28] 

13.  CH4+C+=H+C2Hj 

4.80E+ 14 

0.0 

0.0 

[28] 

14.  CH4+C^=H2+C2Hi 

2.40E+ 14 

0.0 

0.0 

[28] 

15.  C  +  H+=H+C+ 

1.80E+ 15 

0.0 

0.0 

[29] 

16.  0  +  C0+=C0  +  0+ 

8.40E+13 

0.0 

0.0 

[29] 

17.  E+CO+=C+0 

1.86E+18 

-0.5 

0.0 

[29] 

18.  E+CO=E+E  +  CO+ 

5.88E+12 

0.5 

300,000.0 

[15] 

19.  E+CO=E+E  +  0  +  C+ 

1.68E+12 

0.5 

300,000.0 

[15] 

20.  E+CO=E+E  +  C  +  0+ 

8.40E+11 

0.5 

300,000.0 

[15] 

21.  E+C02=E+E+C02 

5.46E+ 12 

0.5 

300,000.0 

[15] 

22.  E+C02=E+E+0+C0+ 

8.40E+11 

0.5 

300,000.0 

[15] 

23.  E+C02=>£+£+02+C+ 

8.40E+ 11 

0.5 

300,000.0 

[15] 

24.  E+C02=E+E+C0+0+ 

8.40E+11 

0.5 

300,000.0 

[15] 

25.  E+C02=E+E+C  +  02 

4.20E+11 

0.5 

300,000.0 

[15] 

26.  CO+E=C+CT 

1.60E+07 

0.0 

0.0 

[15] 

27.  C02+E=CO+CT 

3.00E+11 

0.0 

0.0 

[15] 

28.  0+E=CT 

2.60E+06 

0.0 

0.0 

[15] 

29.  0+E+M=CT+M 

3.60E+16 

0.0 

-600.0 

[15] 

30.  CO + AR+ = AR + CO+ 

4.80E+13 

0.0 

0.0 

[15] 

31.  C++N02  =  C+N02 

1.97E+ 18 

-0.5 

0.0 

[15] 

32.  AR++N02=AR+N02 

5.00E+ 19 

-0.5 

0.0 

[15] 

33.  AR+ + N03  =  AR + N03 

5.10E+17 

-0.5 

0.0 

[15] 

34.  C++NOi=C+N03 

2.00E+ 18 

-0.5 

0.0 

[15] 

35.  C02  +  AR+=AR+C02 

4.56E+ 14 

0.0 

0.0 

[15] 

36.  0+AR+=AR+0+ 

3.84E+ 12 

0.0 

0.0 

[15] 

37.  02+AR+=AR+02 

2.76E+13 

0.0 

0.0 

[15] 

38.  C02+C+=CO  +  CO+ 

6.60E+ 14 

0.0 

0.0 

[15] 

39.  02+C+=C0+0+ 

3.68E+ 14 

0.0 

0.0 

[15] 

40.  02+C+=0  +  C0+ 

2.25E+ 14 

0.0 

0.0 

[15] 

41.  C02  +  CO+=CO+COi 

6.00E+ 14 

0.0 

0.0 

[15] 

42.  02  +  CO+=CO+Oi 

7.20E+13 

0.0 

0.0 

[15] 

43.  0  +  COi=CO  +  Oi 

9.80E+13 

0.0 

0.0 

[15] 

44.  0  +  C0i=C02  +  0+ 

5.77E+13 

0.0 

0.0 

[15] 

45.  02  +  C0i=C02  +  0i 

3.84E+13 

0.0 

0.0 

[15] 

46.  C02  +  0+=C0+0i 

2.70E+ 14 

0.0 

0.0 

[15] 

47.  C02  +  0+=0  +  C0i 

2.70E+ 14 

0.0 

0.0 

[15] 

48.  CO  +  CT=C02+E 

3.30E+14 

0.0 

0.0 

[15] 

constants  measured  at  much  lower  temperatures  (1000- 
3000  K),  appropriate  to  flames.  Some  temperature  dependencies 
of  rate  constants  [6-9]  were  approximated  by  dependencies, 
which  are  valid  at  flame  temperatures.  Their  extrapolation  can 
lead  to  erroneous  results.  The  kinetic  model  was  examined  for 
possible  incorrect  use  of  extrapolated  temperature  dependencies 
of  rate  constants.  Rate  constants  with  large  n  values  (negative  or 
positive)  were  modified  using  literature  data,  where  they  were 
available,  or  were  re-scaled. 

The  Chemkin  package  [3 1  ]  was  used  for  kinetic  modeling  of 
reaction  processes  in  the  plasma  plume.  Reaction  pathway 
analysis  and  calculation  of  sensitivity  coefficients  were  applied 
to  characterize  the  main  reactions  and  the  reaction  mechanism 
in  the  plasma  kernel.  Kinetic  modeling  was  performed  with 
parameters  which  correspond  to  our  experimental  conditions 
(temperature  range,  cooling  rate,  times  of  emission  registration). 
It  was  assumed  that  the  heating  of  the  gas  takes  place 
instantaneously,  since  the  duration  of  the  laser  shot  and  the 
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Table  4 

Time  to  attain  equilibrium  at  a  constant  temperature  (10  atm,  mixture  0.2  O2+O.O3 
H2O+0.77N2) 


Temperature 

Time  (model  without  ionization) 

Time  (model  with  ionization) 

10,000  K 

2  x 10~5 

2  x  10~5 

12,000  K 

7x  10~6 

~7X  10~6 

15,000  K 

2.7  x  10~6 

3.5  x  10~6 

20,000  K 

2.5  x  10  6 

~3  x  10~6 

30,000  K 

2.3  x  10~6 

~3  x  10~6 

generation  of  the  plasma  is  short  compared  with  the  beginning 
and  duration  of  the  LIBS  signal.  Assumed  initial  temperatures 
were  in  the  range  12,000-20,000  K.  These  temperatures 
correspond  to  experimentally  observed  temperatures  [32,33],  It 
was  assumed  that  the  cooling  of  the  gas  follows  exponential 
temperature  decay  with  the  characteristic  decay  times  of  2- 
30  ps,  corresponding  to  experimental  measurements  and  to 
typical  times  of  LIBS  registration  of  the  emission  from  excited 
atoms  in  the  plasma  kernel  [32-34].  Note  that  in  real 
experimental  conditions,  pressure  is  a  variable  parameter  and, 
for  the  studied  time  range,  the  sample  gas  will  possibly  remain 
largely  unmixed  with  a  buffer  gas.  Nevertheless,  even  taking  into 
account  these  more  realistic  conditions,  the  underlying  features 
of  chemistry  and  the  mechanism  of  reaction  process  should 
remain  the  same.  Below  a  short  discussion  will  be  presented 
concerning  the  temperature  dependence  of  reaction  rates, 
reaction  times  and  the  level  of  equilibrium  electron  densities 
for  LIBS  plasma  conditions. 

2.1.  Leveling-off  of  temperature  dependence  of  reaction  rate. 
Reaction  time 

It  should  be  noted  that  temperatures  of  the  order  of  15,000  K 
(and  higher),  appropriate  to  LIBS  plasmas,  actually  lead  to 
relatively  similar  rate  constants  for  different  reactions,  since  the 
Arrhenius  exponential  temperature  dependence  of  reaction  rates 
levels  off  at  high  temperatures.  Thus,  the  characteristic  reaction 
times  for  radicals  and  stable  molecules  become  similar,  even 
with  differing  activation  energies  of  reaction.  Even  with  a  further 
increase  in  temperature,  a  limiting  reaction  time,  e.g.  the  time 


required  to  attain  equilibrium,  is  observed.  Thus  a  high-tem- 
perature  decomposition  of  typical  explosives  very  quickly  leads 
to  a  set  of  simple  species  consisting  of  N,C,H  and  O  atoms.  The 
largest  time  of  decomposition  of  species  with  two  or  three  atoms 
belongs  possibly  to  the  nitrogen  molecule.  We  estimate  that  the 
time  to  attain  equilibrium  in  the  N2  system  is  2-4  ps  at  1 5,000  K 
and  10  atm.  The  time  to  reach  equilibrium  conditions  is  likely  in 
the  ps  range,  and  an  increase  in  temperature  should  not 
significantly  affect  this  time.  Note  that  inverse  bremsstrahlung 
and  multiphoton  ionization  processes  (which  were  not  con¬ 
sidered  in  this  work)  will  accelerate  chemical  reaction  processes 
during  the  initial  period.  However,  the  relatively  slow  reactions 
will  determine  the  time  to  attain  equilibrium. 

Calculated  times  to  attain  thermodynamic  equilibrium  at 
different  temperatures  (under  constant  temperature)  are  pre¬ 
sented  in  Table  4  for  the  N/H/O  system.  The  initial  gas  mixture 
composition  was  0.202  +  0.03H20+0.77N2.  This  composition 
corresponds  to  humid  air  atmosphere  and  includes  gases  which 
are  typically  found  in  the  products  of  decomposition  of  solid 
explosives.  The  results  may  lead  to  several  conclusions.  The 
relatively  long  time  to  attain  equilibrium  suggests  that  the  initial 
decomposition  chemistry  of  explosives  may  be  of  importance  in 
the  chemistry  of  the  LIBS  plume.  Also,  it  means  that  the 
chemistry  may  influence  the  concentration  level  of  emitting 
species,  leading  to  a  chemical  effect  on  the  LIBS  signature.  Of 
course,  the  temperature  and  pressure  conditions,  which  affect 
electronic  quenching,  are  the  most  important  factors  in 
determining  the  concentration  levels  of  the  emitting  species. 
Nevertheless,  as  it  will  be  demonstrated  below,  the  LIBS  plasma 
loses  memory  of  the  initial  compound  at  substantially  shorter 
times  in  comparison  with  equilibration  times. 

In  calculations  the  equilibration  time  was  determined  as  the 
time  to  reach  95%  of  the  equilibrium  concentration.  The  quoted 
times  represent  estimates  under  the  assumption  of  our  kinetic 
model  and  the  rate  constants  employed  in  the  model. 

2.2.  Equilibrium  electron  density  in  different  gases 

The  importance  of  electron-impact  processes  can  be  seen  from 
the  calculated  equilibrium  electron  densities.  Fig.  1  displays  the 


Temperature,  K 
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Fig.  1 .  Temperature  dependence  of  the  equilibrium  electron  density  in  different  buffer  gases  (5  atm).  Calculations  were  performed  with  the  use  of  database  and  program  [30]. 
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temperature  dependence  of  the  equilibrium  electron  density 
calculated  for  several  buffer  gases  and  RDX.  Results  were 
obtained  for  pure  gases.  The  plot  demonstrates  the  electron 
densities,  which  may  develop  within  the  LIBS  plasma  plume. 
Calculations  were  conducted  with  the  use  of  a  database  and 
program  [30],  The  highest  equilibrium  electron  concentration  is 
observed  in  an  atmosphere  of  xenon,  and  the  lowest  equilibrium 
concentration  corresponds  to  helium.  The  differences  in  the 
electron  density  in  different  buffer  gases  are  mainly  determined  by 
differences  in  ionization  potentials.  Results  show  that  temperature 
dependence  of  equilibrium  electron  concentration  for  RDX  gas 
resembles  the  dependence  for  the  electron  density  in  nitrogen.  The 
high  levels  of  electron  densities  demonstrate  that  contributions  of 
electron-impact  processes  are  substantial. 

3.  Results  and  discussion 

Fig.  2  presents  the  temporal  profiles  of  the  concentrations  of 
the  main  species  in  the  RDX  LIBS  plume  for  a  typical  set  of 
initial  conditions.  Due  to  the  high  initial  temperature,  thennal 
reactions  lead  initially  to  the  formation  of  metastable  and 
excited  atoms,  which  are  ionized  in  thermal  processes.  The 
formation  of  electrons  in  thennal  processes  leads  to  a  further 
rapid  accumulation  of  ions  and  electrons  through  electron- 
impact  ionization  processes.  For  species  with  high  ionization 


Time,  s 


Fig.  2.  Time  dependence  of  concentrations  of  the  main  species  in  a  plasma 
plume  of  RDX  (model  of  initial  RDX  decomposition  [8]  with  modifications; 
initial  temperature  18,000  K,  pressure  10  atm,  1%  of  RDX  in  Ar,  characteristic 
temperature  quenching  time  25  p,s). 


potential,  ionization  mainly  occurs  via  electron  impact.  Atomic 
carbon  has  a  relatively  low  ionization  potential.  The  ion  C+  can 
be  produced  by  electron-impact  ionization 

C  +  e— »C+  +  e  +  e  +  266kcal/mol. 

Additionally  charge  transfer  reactions,  which  take  place  in 
the  RDX  LIBS  plume,  contribute  to  C+  formation, 

CO+  +  C->-CO  +  C+-63kcal/mol 

NO+  +  C— >NO  +  C+-46kcal/mol 

Cfy  +  C-»02  +  C+-19kcal/mol 

Fig.  2  contains  time  dependencies  of  the  concentrations 
of  excited  emitting  states  of  oxygen  (777  nm)  and  nitrogen 
(746  nm)  atoms.  Reaction  pathway  analysis  demonstrates  that 
the  main  formation  reactions  of  excited  states  are  electron- 
impact  excitation  from  the  lower  electronic  states.  Electron- 
impact  dissociation  of  02  and  N2  makes  a  small  contribution 
due  to  the  high  degree  of  dissociation.  The  main  processes 
removing  the  excited  states  are  also  electron-induced  quenching 
processes.  Sensitivity  calculations  reveal  that  these  excited 
states  are  also  sensitive  to  the  rate  constants  of  the  following 
reactions  of  neutral  species:  02  +  M=0  +  0  +  M;  OH+M  =  0  + 
H+M;  N2+M=N  +  N  +  M;  N  +  NO=N2  +  0;  NH+H=N+H2; 
CN  +  M=C+N+M;  NCO  +  H=NH  +  CO;  NCO  +  0=NO  +  CO; 
CO  +  M=C+0+M.  The  emission  intensities  can  be  calculated 
by  multiplying  the  concentration  of  emitting  states  by  the 
radiative  transition  probabilities  within  the  framework  of  our 
kinetic  model. 

Kinetic  modeling  of  reactions  in  the  LIBS  plume  demonstrates 
that  the  reaction  mechanism  can  be  considered  to  occur  in  two 
stages.  In  accord  with  our  previous  results  [4,5]  the  first  stage  is 
characterized  by  the  approach  to  the  quasi-stationary  conditions. 
The  high  electron  concentration  sustains  high  populations  of  the 
excited  states  of  H,  N  and  O  atoms.  The  duration  of  this  stage  is 
less  than  ~  1  ps.  During  the  second  stage,  the  temperature  starts  to 
decrease  substantially,  and  the  concentrations  of  species  are  in 
quasi-equilibrium  at  the  decreasing  temperature. 

Analysis  of  the  two-stage  reaction  behavior  shows  that  the 
difference  in  the  mechanisms  of  initial  RDX  decomposition  is 
important  only  during  the  first  reaction  stage.  Fig.  3  contains  the 
results  of  calculations  for  two  cases  of  different  decomposition 
kinetics  of  RDX:  a)  gas-phase  decomposition  of  RDX  in  argon 
using  the  decomposition  kinetic  equation  [10],  was  applied;  and 
b)  gas-phase  product  distribution  at  the  RDX  surface  under  laser 
irradiation,  measured  by  Fetherolf  and  Litzinger  [11]  in  argon, 
was  taken  as  the  initial  conditions.  The  results  after  about  1  ps 
show  very  similar  reaction  behavior.  This  allows  one  to  simplify 
the  kinetic  model  for  calculations  of  reaction  behavior  in  the  LIB  S 
measurement  window  (1-30  ps).  The  results  of  the  analysis  of 
this  two-stage  reaction  behavior  suggest,  that  for  modeling  of  the 
reaction  proceeding  in  the  LIBS  measurement  window  —  1— 
50  ps,  the  knowledge  of  initial  reactions  of  explosive  decom¬ 
position  is  not  important.  As  a  result  of  this  behavior,  the 
developed  model  can  be  used  to  describe  a  reaction  behavior  of 
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Time,  s 

Fig.  3.  Time  dependence  of  O,  N,  N2,  N+,  0*6  (777  nm)  concentrations 
calculated  for  two  models  of  RDX  decomposition:  (1)  decomposition  equation 
from  [10];  (2)  the  set  of  gas-phase  decomposition  products,  measured  at  RDX 
surface  under  laser  irradiation  [11],  was  taken  as  the  initial  mixture  composition 
in  Ai*  (1.38%  RDX  in  Ar,  other  conditions  remain  as  in  Fig.  2). 

other  C/H/N/O  explosives  in  the  gas-phase  LIBS  plume 
environment. 

Fig.  4  contains  results  of  kinetic  calculations  of  the  plasma  with 
four  different  initial  mixtures  corresponding  to  the  same  N/O/H 
atomic  composition.  It  demonstrates  that  a  relatively  fast  approach 
to  approximately  the  same  quasi-stationary  state  is  observed 
between  10  8  and  10~7  s.  The  inclusion  of  carbon  chemistry 
increases  this  time  (Fig.  3).  For  the  large  number  of  species 


considered  in  the  model,  quasi-equilibrium  is  attained  at 
approximately  the  same  times  as  for  the  N-O-H  system. 
However  several  species  experience  a  relatively  late  approach  to 
quasi-equilibrium  at  approximately  1  ps  (e.g.  N2  and  N20). 

The  mechanism  of  attaining  similar  mixture  compositions  for 
compounds  with  the  same  atomic  content  in  the  LIBS  plume  is 
possibly  related  to  a  large  difference  in  the  rate  constants  of  radical 
formation  and  ionization  reactions.  To  some  extent  this  behavior 
is  the  opposite  of  the  two-stage  recombination  behavior  observed 
experimentally:  fast  ion  recombination  with  further  slower  atom/ 
radical  recombination.  Here  fast  decomposition  processes  lead  to 
the  formation  of  the  reaction  mixture  consisting  of  relatively 
simple  atomic  and  diatomic  atom  species  mostly  during 
approximately  the  10-100  ns  period.  Because  this  set  of  species 
is  approximately  the  same  for  initial  mixtures  with  the  same 
atomic  composition,  this  leads  to  a  close  reaction  states  for 
different  mixtures.  Further,  slower  ionization  processes  and 
temperature  evolution  lead  to  changes  in  the  quasi-stationary 
state,  which  is  roughly  the  same  for  mixtures  with  the  same  atomic 
composition  during  the  time  window  of  LIBS  measurements.  It 
means  also  that  quasi-equilibrium  or  quasi-stationary  proceeding 
is  attained  at  later  times  than  times  of  attaining  approximately  the 
same  composition. 

It  is  of  interest  to  note  that  calculations  demonstrate  a  peak  in 
C2  species  concentrations  in  the  interval  between  10  and  400  ns 
depending  on  the  decomposition  rate  and  temperature  decay  time 


Fig.  4.  Results  of  calculations  of  N,  N+,  N2  and  N2  concentrations  for  4  different  mixtures  with  the  same  atomic  composition.  Initial  mixture  compositions  are  shown  on 
the  figure. 
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Time,  s 

Fig.  5.  Double-peak  behavior  of  C2  concentration  (initial  temperature  18,000  K,  10  atm,  1%  RDX  in  argon  atmosphere). 


(Fig.  5).  As  it  was  mentioned  above,  C2  species  were  detected  in 
the  LIBS  plume  of  different  explosives  [2].  The  temperature 
decrease  during  the  second  stage  leads  to  a  further  increase  in  the 
C2  concentration  and  to  the  appearance  of  the  second  maximum. 
The  calculations  show  that  the  initial  temperature  decrease  and  the 
increase  in  quenching  rate  favor  C2  fonnation.  Initially  C2  is 
formed  mainly  in  the  reactions 

CO  +  C^C2  +  0+1 14kcal/mol  (1) 

CN  +  C^C2  +  N  +  36kcal/mol  (2) 

and  in  recombination  processes.  The  slower  dissociation  reactions 
at  high  temperatures, 

CN  +  M— »C  +  N  +  M  +  180kcal/mol 
C2  +  M— >C  +  C  +  M  +  145kcal/mol, 

lead  to  the  appearance  of  the  C2  concentration  maximum.  A 
further  temperature  decrease  leads  to  the  substantial  decrease  in 
the  dissociation  rate,  the  increase  of  recombination  of  C+  ions  and 
the  increase  in  the  concentration  of  C2  that  is  later  consumed  in  the 
reactions  with  radicals. 

Harilal  et  al  [35]  experimentally  observed  two  peaks  in  the  C2 
concentration  during  the  evolution  of  laser  ablated  plasma  from  a 
graphite  target.  They  suggested  that  the  first  maximum  was  the 
result  of  recombination  processes  under  their  experimental 
conditions.  The  second  maximum  was  attributed  to  a  dissociation 
of  carbon  clusters,  leading  to  the  generation  of  C2  molecules  [35], 
Experimental  times  of  the  first  C2  peak  correspond  to  the  times  of 
modeling  results  of  this  work,  which  are  in  the  50-300  ns  range. 
The  second  maximum  was  experimentally  observed  in  the  range 
0.5-3  ]i s.  Modeling  results  show  this  time  is  between  10  and 
50  ps,  depending  on  initial  RDX  concentration  and  characteristic 
temperature  quenching  time.  Of  course,  the  mechanism  of 
appearance  of  the  two  peaks  of  the  C2  species  in  laser  ablation 
of  graphite  may  be  different.  The  appearance  of  second  maximum 
of  C2  species  is  described  by  gas-phase  reactions  in  this  work. 

4.  Conclusions 

In  this  work  we  have  developed  a  kinetic  model  to  describe  the 
processes  which  lead  to  the  electronically  excited  states  of 


oxygen,  hydrogen  and  nitrogen  atoms  responsible  for  the 
emission  lines  in  the  LIBS  plasma  plume  of  RDX  explosive. 
Kinetic  calculations  were  performed  with  the  use  of  temperature 
profiles  with  the  typical  decay  times  in  the  LIBS  plasma.  Reaction 
processes  in  the  RDX  plasma  were  analyzed  through  kinetic 
modeling.  The  main  generation  reactions  of  excited  states  are 
electron-impact  processes.  It  was  found  that  the  evolution  of  the 
concentration  of  C2  species  in  the  RDX  plasma  plume  demon¬ 
strates  double-peak  behavior.  The  mechanism  of  this  bimodal  C2 
concentration  profile  was  discussed.  It  is  of  interest  that  double¬ 
peak  behavior  of  C2  species  was  experimentally  observed  during 
laser  ablation  of  graphite  by  Harilal  et  al.  [35],  Experimental  times 
of  the  observed  C2  peaks  approximately  correspond  to  the 
modeling  results. 

It  was  demonstrated  that  the  overall  reaction  process  in  the 
plasma  plume  can  be  considered  as  a  two-stage  process.  In 
agreement  with  our  previous  results  (Pb/Ar/air  model  [4,5]), 
the  first  stage  corresponds  to  a  fast  approach  to  a  quasi¬ 
stationary  state.  The  second  stage  corresponds  to  a  change  in 
the  quasi-stationary  concentrations  of  species  as  a  result  of  the 
change  in  temperature.  As  a  result  of  the  two-stage  process,  the 
initial  composition  at  the  surface  of  the  target,  or  the  initial 
mechanism  of  RDX  decomposition,  is  not  important  in 
determining  the  LIBS  signature  in  the  LIBS  measurement 
time  window.  The  results  of  modeling  with  different  mechan¬ 
isms  of  decomposition  are  approximately  the  same  during  the 
second  stage.  The  memory  of  the  initial  compound  is  lost 
during  the  period  of  approximately  1  ps.  This  allows  one  to  use 
the  suggested  model  for  modeling  of  the  LIBS  plume  of  other 
C/N/H/O-containing  energetic  materials  in  the  LIBS  measure¬ 
ment  window  (1-30  ps).  Note,  that  the  identification  of 
explosives  is  possible  using  their  unique  ratios  of  atomic 
intensities  [2], 
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